Astronomy & Astrophysics manuscript no. 0388 


© ESO 2008 


August 7, 2008 





On the equilibrium rotation of Earth-like extra-solar planets 

Alexandre CM. Correia 1 ' 2 , Benjamin Levrard 2 ' 3 , and Jacques Laskar 2 



OO 

o 
o 

(N 

< 



Oh 
6 



o 



OO 
O 
OO 

o 



1 Departamento de Ffsica da Universidade de Aveiro, Campus Universitario de Santiago, 3810-193 Aveiro, Portugal 

2 Astronomie et Systemes Dynamiques, IMCCE-CNRS UMR 8028, 77 Avenue Denfert-Rochereau, 75014 Paris, France 

3 Universite de Lyon 1, CRAL, Ecole Normale Superieure de Lyon, 46 allee d'ltalie, 69364 Lyon Cedex 07, France 



Received ?? May, 2008; accepted ? 



99 



ABSTRACT 



The equilibrium rotation of tidally evolved "Earth-like" extra-solar planets is often assumed to be synchronous with their orbital mean 
motion. The same assumption persisted for Mercury and Venus until radar observations revealed their true spin rates. As many of 
these planets follow eccentric orbits and are believed to host dense atmospheres, we expect the equilibrium rotation to differ from the 
synchronous motion. Here we provide a general description of the allowed final equilibrium rotation states of these planets, and apply 
this to already discovered cases in which the mass is lower than 12 M e . At low obliquity and moderate eccentricity, it is shown that 
there are at most four distinct equilibrium possibilities, one of which can be retrograde. Because most presently known "Earth-like" 
planets present eccentric orbits, their equilibrium rotation is unlikely to be synchronous. 
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1. Introduction 

After a significant number of discoveries of new extra-solar 
gaseous giant planets, a new barrier has been passed with the de- 
tections of several planets in the Neptune and even Earth-mass 
(Af ffi ) regime: 5 - 12 M m dRiyera et al.ll2005h lLovis et al.ll2006t 
lUdrv et al.ll2Tj07tlBonnls et al.ll2007l) . If the commonly accepted 
core-accretion model can account for the formation of these 
planets, resulting in a mainly icy/rocky composition, the fraction 
of the residual He-H2 atmospheric envelope accreted during the 
planet migration is not ti ghtly constrained fo r planets more mas- 
sive than the Earth (e.g. lAlibert et al.ll2.Q06h . A minimum mass 
of below 10 M @ is usually consider ed to be t he bo undary be- 
tween terrestrial and giant planets, but Rafikov (2006) found that 
planets more massive than 6 M® could have retained more than 
1 M @ of the He-H2 gaseous envelope. For comparison, masses 
of Earth's and Venus' atmosphere are respectively ~ 10~ 6 and 
10~ 4 times the planet's mass. Despite significant uncertainties, 
these discoveries of "super-Earths" provide an opportunity to 
test some properties that could be similar to those of our more 
familiar telluric planets. 

Becaus e some of these p l anets are potentia lly in the "habit- 
able zone' ' (lUdry et al.ll2007l ISelsis et al.ll2007l) . their spin state 
is an important factor in understanding their climate. For plan- 
ets with a negligible atmosphere such as Mercury, solid tides 
induced by the host star are expected to despin the planet un- 
til an equilibrium position or a capture in a spin-orbit reso- 
nance, depending on the orbital ecc entricity and the perma- 
nent quadrupole mome nt of inertia dGoldreich & PealeHl966t 
ICorreia & Laskar 2004). However, for planets with dense atmo- 
sphere such as Venus, thermal atmospheric tides, driven by solar 
insolation, have a profound influence on the sp in, and may desta- 
bilize the pre vious tidal equilibrium state. ICorreia & Laskatj 
(20011 l2003bl) investigated the combined effect of gravitational 
and atmospheric tides on Venus and showed that the planet could 
evolve into two different final states, the existing retrograde mo- 



tion being the most probable one. This study was performed with 
the small eccentricity approximation that may no longer be ad- 
equate for "Earth-like" planets, which exhibit a wide range of 
eccentricitie s, orbital distances, or ce ntral star types. Following 
the work of lLaskar & Correial (120041) . we generalize here these 
previous studies and investigate the possible equilibrium rotation 
states of "Earth-like" planets in the presence of significant ec- 
centricity. Although our knowledge of these planets is restricted 
to their orbital parameters and minimum masses, we attempt to 
place new constraints on their surface rotation rate, assuming 
that they have a dense atmosphere. 

2. Tidal evolution 

Tidal effects are generated by differential and inelastic defor- 
mation of a planet by a perturbing body. For a planet with a 
dense atmosphere, we consider the traditional gravitational tides 
(t = g), and also thermal atmospheric tides (r = a) produced by 
the stellar heating of the atmosphere. In both cases, the averaged 
variation in the spin can be expressed by the variables that char- 
acterize the spin (the rotation rate a> and the obliquity e), and 
the e lliptical elements, with semi-major ax is a and eccentricity e 
(e.g. lKaulal 19641 ICorreia & Laskdl2003ih : 

dL 

dt 



= r T = J s: T ^M^)A T , £r (e,£), 



(1) 
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where L = Ca> is the angular momentum of the planet, C is the 
polar moment of inertia, the sum being taken over different har- 
monics of the tidal frequency <x (an integer combination of the 
rotation rate a> and the mean motion n), K T is a factor related 
to the strength of the tide, b T {cr) is an odd function related to 
the dissipation within the planet, and A TjCr are polynomial func- 
tions of e and cos e. Since discovered Earth-like planets have 
moderate eccentricities (typically e < 0.25), we neglect terms in 
e 4 , which allows us to simplify the tidal torques acting on the 
planet. In the same way, we consider that the planet's obliquity 
e is small, neglecting terms higher than s 2 in Eq.(Q~|i. 
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Gravitational tides are raised on the planet by the star be- 
cause of the effect of the g ravitational g radient across the planet. 
Equation (fTJ becomes (eg. lKaulall964i) : 



2 8 



(l -5e 2 ) b g (2aj-2n) 



1 

4 



■ 2 49 - 

+ — e b g (2(jj — n) + — e b g (2a> — 3n) 



(2) 



where K g - -GM 2 R 5 /a 6 , G is the gravitational constant, M, is 
the mass of the star, and R is the mean radius of the planet. Since 
planets are not perfectly rigid, there will be a distortion that pro- 
duces a tidal bulge of amplitude k%, the second order potential 
Love number. Imperfect elasticity delays the planet's response 
to the perturbation by a time lag At g . The deformation there- 
fore lag behind the perturbation by an angle 5 g = crAt g /2 and 

b g (cr) = k 2 sm25 g = k 2 sin (cr At g ^ (e.g. lCorreia & Laskarll2003al) . 

The differential absorption of the stellar heat by the plan- 
etary atmosphere gives rise to local variations in temperature 
and consequently to pressure gradients. The mass of the atmo- 
sphere is then redistributed, adjusting for an equilibrium posi- 
tion. Observations on Earth show that the pressure redistribution 
is essentially a superposition of two pressure waves: a diurnal 
tide of small amplitude (wit h no dynamical counterpar t) and a 
strong semi-diurnal tide (see Chapman & Lindzenl 19701) . As for 
gravitational tides, the redistribution of mass in the atmosphere 
produces an atmospheric bulge that modifies the gravitational 
po tential generated by the atmosphere. The resulting tidal to rque 
is dDobrovolskis & IngersoHHl980tlCorreia & Laskail2003al) : 

T a = \k„ [(l - 3e 2 ) b a {2u - 2n) 

- e 2 b a (2co - ri) + 9 e 2 b a (2u> - 3n)] , (3) 

where K a — -3M„R 3 /(5pa 3 ), and p is the mean density of the 
planet. The amplitude of the bul ge is given by p?, the seco nd 
order surface pressure variations (Chap man & Lindzenlll970l) : 



p 2 (cr) =i-p Q \V-V (r 

cr 



y-i Jo- 



y gH ) cr 



.Pa 



(4) 



where y = 7/5 for a perfect diatomic gas, pa is the mean surface 
pressure, v is the velocity of tidal winds, J is the amount of heat 
absorbed or emitted by a unit mass of air per unit time, and Hq 
is the scale height at the surface. There is also a delay At a be- 
fore the response of the atmosphere to the stellar heat excitation. 
However, the imaginary number in Eq.© causes the pressure 
variations to lead the Sun (i = e"' 2 ). Using 5 a - crAt a /2, we thus 
have b a (cr) = p 2 sin2£ fl = \p 2 \ s\r\2{5 a + n/2) = -\p 2 \ sin2tf a . 

The dependence of time lags on the tidal frequency is poorly 
known for gravitational and atmosph eric tides. For a revi ew of 
the different dissipation models, see ICorreia et al.l d2003l) . For 
both tides we adopt the simplest mo del for slow rot ation (i.e. 
the viscous model), as described by Mignardl d 19791) . where a 
constant time lag is assumed for all components of the tidal 
perturbation. Since we usually have crAt <k 1, this model can 
be made linear: b g (cr) - kicrAt g and b a (cr) - -\p2\ crAt a . 
For atmospheric tides, it is also necessary to consider the re- 
sponse of the surface pressure variations to tidal frequency 
(Eq.|4]i. We use the "heating-at-t he-ground model" described by 
Dobrovolskis & IngersolH( ll980h . It is supposed that all the stel- 
lar flux absorbed by the ground, F s , is immediately deposited 
in a thin layer of atmosphere at the surface. The heating dis- 
tributing is then written as a delta-function just above the ground 



(Jo- - gFJpo). This approximation is justified because tides in 
the upper atmosphere are decoupled from the ground by the dis- 
parity between their rotation rates and it appears to be in good 
agreement with the observations. Neglecting v a over the thin 
heated layer, Eq.(@) becomes: Va- = F s /(8Ha) oc L*/a 2 (L„ is 
the star luminosity). 

The average evolution of the rotation rate is obtained by 
adding the effects of both tidal torques acting on the planet, that 
is lo — (T a + T g )/C. Substituting dissipative models described 
above in Eqs.(f2| and (f3]) then leads to 

o)/Kq = a) — ( 1 + 6<? 2 )n - a> s [(l - 3e 2 ) sign(w - n) 

- e 2 sign(2w - ri) + 9 e 2 sign(2w - 3n)] , (5) 

where K = -m(e)K g k 2 At g /C, Cl(e) = 1 + 15e 2 /2, 



K a At a 



R 

— a 



8H £l(e)k 2 K g At„ AT, m 



(6) 



and m is the planet mass. According to Eq.©, atmospheric 
tides are weak during the first stages of evolution (cr » ri). 
Gravitational tides alone can therefore be used to estimate the 
characteristic time needed to reach the equilibrium rotation, T eq : 



3k 2 At g R 5 9GMlk 2 At g R 3 



CG 



(7) 



with C - mR 2 /3. All terrestrial extra-solar planets listed in 
Table[T]have a despinning timescale that is significantly lower 
than the age of the system (< 10 7 yrs), so we expect that they 
have already reached their equilibrium rotation state. 



3. Equilibrium final states for the rotation rate 

An equilibrium final state is achieved when co — 0, that is for 
T„ = —T a . When e — 0, we derive: 



f(io -«) = — 



K g b g (2u> — 2ri) 



1 



(8) 



T a K a b a (2(L> - 2ri) 

where f(x) is an even function of x (ICorreia &Laskarfl200lh . 
Assuming that fix) is monotonic close to the equilibrium (which 
is true for the usual dissipation models), we have 



\u-n\=r\\) = a> s , 



(9) 



i.e. there are two final possibilities for the equilibrium rotation of 
the planet, given by cj* = n±o) s , where 2n/cj s can be seen as the 
synodic period. Unless a> s <k n, the synchronous motion (a> - ri) 
is no longer the final equilibri um position for Earth-like planet s 
tidally evolved. According to ICorreia & Laskarl (1200 lL [2003b), 
the planet is free to evolve to any of the two final equilibrium 
positions, although it appears to have a preference for the slow- 
est one {of = n — cj s ). If u> s < n, both final states correspond to 
prograde final rotation rates (FigQJ),c). However, if a> s > n, the 
slowest rotation rate becomes negative, corresponding to a ret- 
rograde configuration. This is the present situation of the planet 
Venus for which o) s /n = 1.92 (FigQJ). 

When e + 0, Eq.© is no longer valid and additional equi- 
librium positions for the rotation rate may occur. For moderate 
values of the eccentricity, we can write from Eqs.(f2]i and (0: 



f(co-n) = 1 +e z 2 + 



b a (2iL) - 2ri) 



(10) 
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Fig.l. Evolution of 0) (Eq.EJl with (o s /n = 1.92, cojn = 0.40 
and cojn = 0.05, for different eccentricities (e = 0.0,0.1,0.2). 
The equilibrium rotation rates are given by u> = and the arrows 
indicate whether it is a stable or unstable equilibrium position. 
For a> s /n > 1, we have two equilibrium possibilities, o>J, one of 
which corresponds to a retrograde rotation (as for Venus). For 
tOs/n < 1, retrograde states are not possible, but we can still ob- 
serve final rotation rates aT < n. For eccentric orbits, because of 
the terms in b T {2io - n) and b T (2a> - 3n), we may have at most 
four different final possibilities (EqfTZl). When a> s /n becomes ex- 
tremely small, which is the case for the present observed extra- 
solar planets with some eccentricity (TableQ]), a single final equi- 
librium is possible for a»|. 



where 



K 



g(<o) = j§r [b g (2co-n)+ 49b g (2<o - 3n)] 
- b a (2u) - n) + 9 b a (2oj - 3n) . 



(id 



Since b T (cr) are monotonic odd functions, the effect of the eccen- 
tricity is eventually to split each previous equilibrium rotation 
rate into two new equilibrium values so that four final equilib- 
rium positions for the rotation rate are possible, written as: 



a>t 9 = n ± oj s + e 2 5* 



'1,2 " 
with 



1,2 ' 



6* - ( 2 1 gM \ ° r> 
u \ \b a {2co-2n)\) dx 

or, adopting the tidal models described in Sect. |2](Eq.[5]): 

8~ 2 = 6n - (6 + 1) w, and 5\ 2 = 6n - (4 + 9) a> s , 



(12) 



(13) 



(14) 



where + corresponds to the state 5\ and - to the state 62- Because 
the set of a>? 2 values must verify the additional condition 



< n/2 < o) l < n < oj\ < 3n/2 < aS^ , 



(15) 



these four equilibrium rotation states cannot, in general, exist 
simultaneously, depending on the values of u> s and e. In par- 
can never coexist with aC. 



ticular, the final states 



and cjt 



At most three different equilibrium states are therefore possible, 
obtained when a> s /n is close to 1/2, or more precisely, when 
1/2- 17e 2 /2 < w. s /n < 1/2 + 7 e 1 12. Conversely, we found that 
one single final state w[ = ( 1 + 6e 2 ) n + ( 1 - 1 3e 2 ) oj s exists when 
cojn < 6e 2 (l - 7e 2 ). 



4. Application to Earth-like extra-solar planets 

The Earth and Venus are the only Earth-like planets for which 
the atmosphere and spin are known. Only Venus is tidally 
evolved and therefore suitable for applying the above expres- 
sions for tidal equilibrium. We can nevertheless investigate the 
final equilibrium rotation states of the already detected "super- 
Earths". For that purpose, we considered only the 6 extra- 
solar planets of masses smaller than 12M e that we classified 
as rocky planets with a dense atmosphere, although we stress 
that this mass boundary is quite arbitrar ily. Using the emp irical 
mass-luminosity relation L* oc Ml (eg. ICester et af]ll983h and 
the mass-radius r elationship for terrestrial planets R oc m 274 
dSotin et al.ll2007l) . Eq.© can be written as: 



tojn 



k(aM*) 2 - 5 m-°- 726 



(16) 



where A: is a proportionality coefficient that contains all the con- 
stant parameters, but also the parameters that we are unable to 
constrain such as Ho, k%, At g or At a . In this context, as a first or- 
der approximation, we consider that for all these terrestrial plan- 
ets, the parameter k has the same value as for Venus. Assuming 
that the rotation of Venus is pres ently stabilized i n the aT final 
state, that is, 2nj(jjT = -243 days dCarpenterl 19701) . we compute 
2n/aj s = 116.7 days. Replacing it in Eq.dToTl. we find for Venus 
that k = 3.32 m°- 726 M 2 - 5 AU -2 - 5 . We can then estimate the ratio 
a) s /n for all considered extra-solar planets in order to derive their 
respective equilibrium rotation rates (Table[T]i. The number and 
values of the equilibrium rotation states are plotted as a function 
of aM„ for different eccentricities in Fig. [2] All eccentric planets 
have a ratio u jn that is lower than 7 x 10~ 3 , which verifies the 
condition ojjn < 6e 2 (l-7e 2 ). As a consequence, only one single 
final state a>+/n (1 + 6e 2 ) exists, corresponding to the equi- 
librium rotation induced by gravitational tides only (Eq.|2}. The 
main reason is that the effect of atmospheric tides is clearly dis- 
favored relative to the effect of gravitational tides on Earth-like 
planets discovered orbiting M-dwarf stars: The smaller orbital 
distance strengthens the effect of gravitational tides, which are 
proportional to I /a 6 , while the effect of thermal tides varies as 
1 /a 5 . Moreover, the smaller mass of the central star also strongly 
affects the luminosity received by the planet and hence the size 
of the atmospheric bulge driven by thermal contrasts. 

The only exception is GJ 876 d, because its eccentricity is be- 
lieved to be zero (yet to be confirmed). However, the two equi- 
librium rotation states a>* are so close to the mean motion n, that 
the quadrupole moment of inertia (not included in our analysis) 
will probably capture the rotation in synchronous resonance. 
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aM* (All M j aM* (All M ) aM* (AU M ) 

Fig. 2. Equilibrium positions of the rotation rate as a function of the product aM, for three different values of the eccentricity 
(e = 0.0, 0.1, 0.2). Each curve corresponds to a different final state (dotted lines for a>^ and solid lines for Wj). For e (case of 
Venus), we always count two final states that are symmetrical about n. For small values of aM*, the two equilibrium possibilities 
are so close to n that the most likely scenario for the planet is to be captured in the synchronous resonance (case of GJ876d). As 
we increase the eccentricity, we can count at most three final equilibrium rotations, depending on the value of u> s /n (computed from 
Eq.[T6li. When e 0.2, only one equilibrium state exists for aM* < 0.3, resulting from a> s /n < 6e 2 (l - le 2 ). This is the present 
situation of GJ581 d and most of the "Earth-like" extra-solar planets listed in TableQ] 



Table 1. Characteristics and equilibrium rotation rates of Earth-like planets with masses lower than 12 M m (see text for notations). 



Name 


M, 


Age 




m sin i 


a 


e 


ojjn 


2n/n 


1nju) x 


27T/lL> 2 


2tt/w[ 






[M ] 


[Gyr] 


[Gyr] 




[AU] 






[day] 


[day] 


[day] 


[day] 


[day] 


Venus 


1.00 


4.5 


2.3 


0.82 


0.723 


0.007 


1.92 


224.7 




-243 




76.8 


GJ581C 1 


0.31 


4.3 


io- 5 


5.0 


0.073 


0.16 


0.0002 


12.93 






11.2 




GJ 876 d 2 


0.32 


9.9 


io- 8 


5.7 


0.021 





10" 5 


1.9378 


1.9379 




1.9377 




GJ581d' 


0.31 


4.3 


0.04 


7.7 


0.253 


0.2 


0.0026 


83.6 






67.4 




HD 69830 b 3 


0.86 


4-10 


io- 5 


10.2 


0.079 


0.10 


0.0009 


8.667 






8.17 




GJ674b 4 


0.35 


0.1-1 


io- 7 


11.7 


0.039 


0.2 


IO 5 


4.693 






3.79 




HD 69830 c 3 


0.86 


4-10 


io- 3 


11.8 


0.186 


0.13 


0.0069 


31.56 






28.5 





U sing Eq.Q wi t h k 2 = 1/3 and At g = 640s (Earth's values); Referen ces: UllUdrv et al.l i2007l) : m lRivera et ai]d2005l) : mlLovis et al.l d200d) : 
r4l lBonfilsetaLH2007h . Earth-like planets OGLE-2005-BLG-390Lb dBeaulieu et alj2006n and MOA-2007-BLG-192-Lb dBennett et al]|2008h 
have not been included because their despinning timescales r eq ~ 10 2 Gyr are much larger than the age of the Universe. 



5. Discussion and conclusion 

We have derived a simple model that permits the determina- 
tion of the final equilibrium rotation of "Earth-like" planets. Our 
model contains some uncertain parameters related to the dissipa- 
tion within the planets, but we were able to gather all this infor- 
mation in a single parameter, co s . By varying this, we can then 
cover all possibilities for the rotation of terrestrial planets. We 
demonstrated that for a planet of moderate eccentricity and low 
obliquity, at most four final equilibrium positions are possible. 
For eccentricities higher than e ~ 0.5, terms of higher degree 
in e should be considered in Eq.(f2]l that may generate additional 
equilibrium positions. 

Based on the present rotation of Venus, we provided an es- 
timate of oj s for different environments (Eq.[T6l). An important 
consequence is that the ratio a> s /n increases rapidly with the 
semi-major axis and mass of the star because a> s /n oc (aM,) 2,5 . 
The effect of the atmosphere on the equilibrium rotation is there- 
fore more relevant for planets that orbit Sun-like stars at not 
close distances. Unfortunately, this situation is precisely the op- 
posite of that for the discovered Earth-like planets (TableQ]), 
since the radial velocity technique is more sensitive to the de- 
tection of short-period planets. As a consequence, unlike Venus, 
none of these planets can be stabilized with a rotation rate 



oj < n, because for e > 0.1 these final states only exist if 
aM* > 0.2AUM Q (Fig©. The effect of atmospheric tides is 
extremely small (cojri ~ 0) and the equilibrium rotation is es- 
sentially driven by the tidal gravitational torque. However, their 
significant eccentricity (0.1 < e < 0.2) moves the equilibrium 
position away from the synchronous motion. Capture in this res- 
onance is unlikely, but capture in a higher order spin-orbit reso- 
nance cannot be rejected, in particular for GJ581 d and GJ674b, 
depending on the asymmetry of their mass distribution. 
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